D( 18 F,pa) 15 N reaction applied to nova 7-ray emission 
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The 18 F(p,q) 15 reaction is recognized to be one of the most important reactions for nova gamma- 
ray astronomy as it governs the early E < 511 keV gamma emission. However in the nova tempera- 
ture regime, its rate remains largely uncertain due to unknown low-energy resonance strengths. We 
report here the measurement of the D( 18 F,p) 19 F(a) 15 N one-nucleon transfer reaction, induced by a 
14 MeV 18 F radioactive beam impinging on a CD2 target; outgoing protons and 15 N (or a-particles) 
were detected in coincidence in two silicon strip detectors. A DWBA analysis of the data resulted 
in new limits to the contribution of low-energy resonances to the rate of the 18 F(p,cv) 15 reaction. 



PACS numbers: 25.60.Je, 21.10.Jx, 26.30.+k, 27.20.+n 
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Gamma-ray emission from classical novae is domi- 
nated, during the first hours, by positron annihilation 
resulting from the beta decay of radioactive nuclei. The 
main contribution comes from the decay of 18 F (half-life 
of 110 min) and hence is directly related to 18 F forma- 
tion during the outburst Q, 0, 0- A good knowledge of 
the nuclear reaction rates of production and destruction 
of 18 F is required to calculate the amount of 18 F syn- 
thesized in novae and the resulting gamma-ray emission. 
The rate relevant for the main mode of 18 F destruction 
(i.e, through 18 F(p,cv) 15 0) has been the object of many 
experiments 0,0] (see also refs. in 0). However, this rate 
remains poorly known at novae temperatures (lower than 
3.5 x 10 8 K) due to the scarcity of spectroscopic informa- 
tion for levels near the proton threshold in the compound 
nucleus 19 Ne. This uncertainty is directly related to the 
unknown proton widths (T p ) of the first three levels (E x , 
r = 6.419 MeV, 3/2+; 6.437 MeV, 1/2" and 6.449 MeV, 
3/2+) . The tails of the corresponding resonances (at re- 
spectively En = 8, 26 and 38 keV) can dominate the 
astrophysical S— factor in the relevant energy range 0. 
As a consequence of these nuclear uncertainties, the 18 F 
production in novae and the early gamma-ray emission 
used to be uncertain by a factor of «300 3]. Unfortu- 
nately, a direct measurement of the relevant resonance 
strengths is impossible due to the very low Coulomb bar- 
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rier penetrability. Hence, we used an indirect method 
aiming at the determination of the one nucleon spectro- 
scopic factors (S) in the analog levels of the mirror nu- 
cleus ( 19 F) by a neutron transfer reaction: D( 18 F,p) 19 F. 
Assuming the equality of one nucleon spectroscopic fac- 
tors in analog levels, it is possible to calculate the proton 
widths through the relation r p = S x T s p where T s p 
is the single particle width. In this rapid communication 
we present the experimental technique, the results of a 
DWBA analysis and the implications on the 18 F(p,a) 15 
thermonuclear reaction rate. 

The experiment has been carried out at the Centre 
de Recherches du Cyclotron at Louvain-la-Neuve (Bel- 
gium) where we used a 18 F radioactive beam. The 18 F 
was produced through the 18 0(p,n) 18 F reaction, chemi- 
cally extracted to form CH^F molecules, transferred to 
the cyclotron source, ionized to the 2 + state and accel- 
erated to 14 MeV0. A total of 15 bunches of <1 Ci 
of 18 F were produced providing each a mean beam in- 
tensity of 2.2xl0 6 particles per second over a period of 
w2 hours. The beam contamination from its stable iso- 
bar was very small, 18 / 18 F < 10~ 3 . The 18 F beam 
was directed on deuteriated polyethylene (CD 2 ) targets 
of wlOO /ig/cm 2 thickness made at Louvain-La-Neuve. 
For the energy considered here (1.4 MeV in the center 
of mass) , the deuteron and the outgoing proton are both 
below the Coulomb barrier. In a similar transfer reaction 
19 F(d,p) 20 F 7] studied at the same cm. energy, strong 
direct contributions have been observed. Hence, with 
a stripping integral dominated by the best known por- 
tion of the wave-function, accurate extraction of spec- 
troscopic factors is expected. 

The experimental setup is depicted in Figure ^ It 
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consists of two silicon multistrip detectors composed of 
sectors with 16 concentric strips The first detector, 
LAMP, is positioned 9 cm upstream from the target; it 
consists of 6 sectors forming a conical shape covering 
laboratory angles between 115° and 160° (i.e. forward 
center-of-mass angles between 12° and 40°) providing 
a good acceptance for protons in the domain of inter- 
est for the differential cross section. The other detector, 
LEDA, is made up of 8 sectors forming a disk positioned 
40 cm downstream from the target and is used for back- 
ground reduction and normalization. The levels of inter- 
ests are situated high above the alpha emission threshold 
(at 4.013 MeV) and their almost exclusive decay mode 
is through 19 F* — > 15 N+a. Hence, to reduce background, 
we required coincidences between a proton in LAMP and 
a 15 N (or an a-particle discriminated by time of flight) 
in LEDA. Following Monte-Carlo simulations, the pro- 
ton detection efficiency is found to be 27% and is only 
slightly reduced to 19% when the coincidence condition 
is applied. 




FIG. 1: Experimental setup. 

Thanks to the kinematics, only protons and a from 
D( 18 F,p) 19 F and D( 18 F,a) 16 can reach LAMP. Using 
energy, angle and time-of-flight information, the events 
of interest with one proton in LAMP are selected. The 
excitation energy of the decaying 19 F levels can be kine- 
matically reconstructed from the energies and angles of 
the detected protons and the known beam energy. The 
energy calibration of the silicon strips has been done with 
a triple-line alpha source and a small correction (« 1 %) 
on the excitation energy was introduced using a linear fit 
to the well-known energies of some lines in D( 18 F,p) 19 F. 
This correction is needed because the reconstructed en- 
ergy is very sensitive to the precise detector position (a 
1 mm error induces a m 10 keV energy shift). The sta- 
tistical error from the fitting procedure induced by this 
correction is about 5 keV at 6.5 MeV. The corresponding 
coincidence spectrum is shown in Figure where verti- 
cal lines correspond to the known position of the 19 F 



levels |9j populated with low transferred angular momen- 
tum (l < 2). The resolution (w 100 keV FWHM) is not 
sufficient to separate the various levels, but the two 3/2 + 
levels of interest at 6.497 and 6.528 MeV (the analogs 
of the 3/2+ levels in 19 Ne) are well separated from the 
other groups of levels. 

The angular distribution shown in Figure|3is the result 
of a selection made on the 6.5 MeV peak in the coinci- 
dence spectrum. The coincidence efficiency is determined 
with a Monte-Carlo simulation taking an isotropic angu- 
lar distribution for the a-decay of the 19 F. Coincidences 
were initially designed to eliminate 18 0(d,p) 19 events 
induced by a possible isobaric beam contamination. This 
was not essential for this purpose because of the high 
beam purity, but it was useful to reduce the background 
from electronic noise affecting the low proton energy re- 
gion. We checked that the ratio between coincidence and 
single events was well reproduced by Monte-Carlo sim- 
ulations. Rutherford elastic scattering of 18 F on carbon 
from the target, detected in LEDA, provides the (target 
thickness) x (beam intensity) normalization. At such low 
beam intensity, the target stoichiometry is not expected 
to be modified during the experiment. This was con- 
firmed by proton elastic scattering analysis of the targets 
done after the experiment at the Orsay ARAMIS facil- 
ity. The solid lines in Figure correspond to theoretical 
DWBA calculations made with the code FRESCO [hJ in 
the zero-range approximation for different transfered an- 
gular momentum (I = 0, 2). The nuclear potentials were 
taken from Ref. [7j in which a similar neutron transfer 
reaction, 19 F(d,p) 20 F, was studied at the same center- 
of-mass energy (subcoulomb transfer). A typical direct 
mechanism pattern was observed for 1 1 angular distribu- 
tions of 20 F excited states. The comparison between our 
data and the calculations indicates a strong dominance of 
I = transfer for the sum of the contribution of the two 
3/2+ levels and the value obtained for the total spectro- 
scopic factor is Sj_ + S' 2 = 0.21 0|. Using an alternate set 
of optical parameters from a compilation |ll| led to the 
same spectroscopic factor value showing the stability of 
the result with respect to variation of the optical param- 
eters. Although no peak is seen in Figure |21 for the l/2~ 
level (E x = 6.429 MeV, T = 280 keV) due to its large total 
width, it is possible to derive an upper limit for the spec- 
troscopic factor of £3 0.15 assuming an / = 1 transfer 
(see Figure El . The important result of this analysis is 
that the contribution of these two 3/2 + resonances to the 
destruction rate of 18 F cannot be neglected. Sources of 
uncertainty on spectroscopic factors are statistics (negli- 
gible), normalization to elastic scattering data (« 10 %) 
but the major one is the DWBA method itself. This last 
one is difficult to evaluate but using statistics of ratios be- 
tween the directly measured proton widths and the ones 
extracted from a DWBA analysis > we estimate this er- 
ror to 25% (r.m.s.) for spectroscopic factors greater than 
0.01. Larger deviation occurs for smaller spectroscopic 
factor due to the contribution of other mechanisms, but 
the high value of 0.21 provides a safe margin. 
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FIG. 2: Example of a reconstructed 19 F coincidence spectrum showing the two 3/2+ levels of astrophysical interest around 
6.5 MeV of excitation energy. The bold line shows the calculated contribution of the l/2~ level at 6429 keV for S3 = 0.15 (see 
text). This spectrum is limited at low energy by the coincidence condition (no 19 F breakup) and at high energy (Ex ~ 7.3 MeV) 
by the proton energy below electronic threshold. 



e 

1 













1 = S = 0.21 


















1 = 2 S = 0.21 







(x10) 



5 10 15 20 25 30 35 40 45 50 



©cm. (cleg) 

FIG. 3: Comparison between the experimental angular distri- 
bution and DWBA calculations for different transferred an- 
gular momentum. The vertical error bars are only statistical 
whereas the horizontal ones are the angular width of each 
strip as seen from the target. 



To show the improvement brought by this experiment, 
we present in Figure 0] the rates calculated in the same 



way as in Ref. 3j except that we use new spectroscopic 
factors obtained from this experiment and the new value 
for the 330 keV resonance strength 0(^7 = 1.48 ± 0.46 
eV). Assuming the equality of spectroscopic factors be- 
tween analog levels one obtains Si + S2 = S[ + S' 2 = 0.21 
and S3 = S 3 < 0.15 where S[, S' 2 and S 3 correspond 
to the E x , r = 6.497 MeV, 3/2+; 6.528 MeV, 3/2+; 
6.429 MeV, l/2~ levels in 19 F while Si, S 2 and S 3 cor- 
respond to the E x , ,V = 6.419 MeV, 3/2+; 6.449 MeV, 
3/2+; 6.437 MeV, 1/2" analog levels in 19 Ne, respec- 
tively. Taking into account the moderate uncertainty on 
the energy calibration, the data favors a dominant con- 
tribution from the 6.528 MeV level i.e. S[ ~ 0. and 
S' 2 ~ 0.2. However the assignment of the two analog 
levels in 19 Ne is not firmly established. We cannot ex- 
clude inversion or mixing of these two 3/2+ levels. Hence 
to be conservative, for the calculation of the rate, we al- 
low the individual spectroscopic factors to take any value 
compatible with the S\ + S2 = 0.21 and S3 < 0.15 con- 
straints. To show the improvement provided by this ex- 
periment, we calculate upper and lower rate limits by 
setting Si = 0, S2 = 0.21 (no inversion), S3 = 0.15 and 
Si = 0.21, S2 = (inversion), S3 = respectively. Here 
we neglect the error associated with the spectroscopic 
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FIG. 4: Present limits to the 18 F(p,a) 15 reaction rate, 
calculated from this work and the new 330 keV reso- 
nance strength^] (hatched area), compared to the previous 
situation 3] (CHJT00). All rates are normalized to the Coc 
et al. nominal rate[IJ used in the more recent nova calcula- 
tions. The Wiescher and Kettner rateO (WK82) used in the 
earliest nova calculations 0] is also displayed for comparison. 



factor extraction discussed above. As a consequence, the 
rate uncertainty is reduced by a factor of w 5 (Figure 
in the temperature range of novae. This is mainly due 
to the reduced contribution of the 1/2" level while the 
impact on the rate of the new contribution of the two 
3/2 + levels is more important at low temperatures. It is 
worth noticing that these new rate limits still encompass 
the nominal rate 3] previously used in the most recent 
gamma-ray flux calculations. Since the repartition of 
spectroscopic factors between the two 3/2 + levels is not 
settled, the choice S\ — 0.1 and 5*2 = 0.1 for the nominal 
rate 3] is still acceptable. However, the new strength^ 
of the resonance at 330 keV has to be used for the rate 
calculation. 
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